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ABSTRACT 

A time series of soft X-ray emission observed on 1974-2007 years (GOES) is 
analyzed. We show that in the periods of high solar activity 1977-1981, 1988- 
1992, 1999-2003 the energy statistics of soft X-ray solar flares for class M and C 
is well described by a FARIMA time series with Pareto innovations. The model 
is characterized by two effects. One of them is a long-range dependence (long- 
term memory), and another corresponds to heavy-tailed distributions. Their 
parameters are statistically stable enough during the periods. However, when 
the solar activity tends to minimum, they change essentially. We discuss possible 
causes of this evolution and suggest a statistical model for predicting the flare 
energy statistics. 

Subject headings: Sun: activity — Sun: flares — Sun: X-rays, gamma rays — 
methods: data analysis — methods: statistical 
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1. Introduction 

Individual solar cycles are different in form, amplitude and length. At present, the 
accurate solar data is only available for the most recent three cycles. Understanding the 
long-term solar variability and predicting the solar activity is an actual problem for so- 
lar physics. It is very important to predict the time and strength of such events because 
these disturbances can pose serious threats to man-made spacecrafts, can disrupt elec tronic 



communication channels and can even set up huge electrical currents in power grids (j Clark 



20061 ). It is enough to remind about serious problems with GOES, Deuthsche Telecom, 
Telstars 401, etc. Satellite operators would be glad to escape the unhappy surprise, and 
mission planners are compelled to take into account the future space weather forecast. Not 
only NASA satellites malfunctioned because of the disturbances, but the global positioning 
system was impaired. The cost to the airline industry arose as planes were re-routed to lower 
altitudes, burning more fuel in force of atmospheric drag. 

As the geological records show, the Earth's climate has always been changing. The 
reasons for such changes, however, have always been subject to continuous discussions and 
are still not well understood. In addition to natural climate changes the risk of human 
influence on climate is seriously considered too. An y factor that alters the radiation received 



from the Sun or lost to Space will affect climate. So. lMann et al. I(jl998l ) have clearly detected 
a significant correlation between solar irradiance and reconstructed Northern Hemisphere 
temperature. The statistics indicates that during "Maunder Minimum" of solar activity the 
climate was especially cold, but when the intensity of solar radiance again increased from 
early nineteenth century through to the mid-twentieth century, the period coincides with the 
general warming. This, however, would either imply unrealistically large variations in total 
solar irradiance or a higher climate sensitivity to radiative forcing than normally accepted. 
Therefore, other mechanisms have to be invoked. The most promising candidate is a change 
in cloud formation because clouds have a very strong impact on the radiation balance and 
because only little energy is needed to change the cloud formation process. According to 
satellite records t aken from 1979 to 1992, Earth was 3 % cloudier during solar minima than 



at solar maxima ( Svensmark fc Friis-Christensen Il997h . One of the ways to influence cloud 



formation mi ght be through the cosm ic ray flux that is strongly modulated by the varying 



solar activity. IScafetta fc West 1(120031 ) argue that Earth's short-term temperature anomalies 



inherit a Levy-walk memory component from the intermittence of solar flares. 

The aim of this paper is to present a statistical model for predicting soft X-ray solar burst 
activity in the period of solar cycle maxima. The paper is organized as follows. The random 
features of solar activity is outlined in Section [2l The data set is described in Section [31 In 
Section H] we present the essence of our statistical investigation of SXR flares. Interrelations 
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among statistical flare parameters, such as long-range dependence index, Hurst exponent 
for X-ray flux and their evolution during solar cycles, are analyzed. In the period of strong 
solar activity the index of self-similarity is nearly constant (Section [S]). This feature can be 
used for predicting the power of soft X-ray emission for the 24-th solar cycle near its solar 
activity maximum (2010-2014). The corresponding model is constructed in Section [61 Next, 
we give a summary and discussion of the main results. Finally, the conclusions are drawn in 
Section [71 



2. Randomness in Solar Activity 



The solar 11-year cycle is driven by Sun's magnetic field. The Sun's magnetic field is 
produced by a hydromagnetic dynamo process underneath th e solar surface and is cyclic in 
nature. This fundamental theoretical idea was established by [Parker I(ll955l ). However, only 
within the last few years, theoretical models of the solar dynamo have become sophisticated 
enough to explain various aspects of the solar activity. So, recently [Pikrati fc Gilman I(l2006l ) 
have made the first attempt of using a theoretical dynamo model to predict the strength 
of the upcomin g cycle 24. They have shown that this cycle will be the strongest in 50 
years. But later [Choudhuri et al. [(|2007[ ) pointed out that some assumptions in the Dikpati- 
Gilman model are uriiustifi ed. On the contrary, their model, based on the earlier work of 
Nandy fc Choudhuri [(|2002[ ). predicts that the cycle 24 will be weaker than the 23-rd. The key 
problem here is the following: the dominant processes like the magnetic field advection and 
toroidal field generation by differential rotation are fairly regular during the rising phase of a 
cycle from a minimum to a maximum, and hence a good knowledge of magnetic configurations 
during a minimum would enable a good theoretical model to predict the next maximum 
reliably. However, the dominant process in the declining phase of a cycle contains the poloidal 
field generation by the Babcock-Leighton mechanism which involves randomness (primary 
cause of solar cycle fluctuations) and cannot be predicted in advance by any deterministic 
model. That is why, although active regions appear in a latitude belt at a certain phase 
of the solar cycle, where exactly within this belt the active regions appear seems random. 
Since the poloidal field generated from an active region depends on the tilt, the scatter in 
the tilts introduces randomness in the poloidal field generation process. 

The other feature of solar activity is that there is a "magnetic persistence" between 
the surfa ce polar fields and sp ot-producing toroidal fields, generated by differential rotation 
shearing (iDikpati et al. This means that the Sun retains a memory of its magnetic 

field for a long time (about 20 years or so). The solar cycle prediction is similar to that 
employed in global atmospheric dynamics over the last ten years. Such models predict 
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changes in certain global characteristics of a cycle, without attempting to reproduce details 
that occur on smaller spatial scales and shorter time scales. The interrelation between global 
characteristics and small scale processes is an open problem, and meanwhile some effects of 
smaller scales are included in parametric form. 



3. X-ray Flare Observations 

Solar activity is a many-sided phenomenon. It includes flares, prominence eruptions, 
coronal mass ejections, solar energetic particles, various radio bursts, high-speed solar wind 
streaming from coronal holes, etc. Solar flares are the most energetic and violent events 
occurring in the solar atmosphere. The energy release in a flare ranges from 10^6 to 3x10^2 
ergs. Magnetic reconnection is considered to play a central role in any flare energy release. 

Observations of solar flare phenomena in X-rays became possible in the 1960s with 
the availability of space-borne instrumentation. Since 1974 broad-band soft X-ray emission 
of the Sun has been measured almost continuously by the meteorology satellites operated 
by NOAA so as the Synchronous Meteorological Satellite (SMS) and the Geostationary 
Operational Environment Satellite (GOES). The first GOES was launched by NASA in 
1975, and the GOES series extends to the currently operational GOES 11 and GOES 12. 
From 1974 to 1986 the soft X-ray records are obtained by at least one GOES-type satellite; 
starting with 1983, data from two and even three co-operating GOES are generally available. 
The X-ray sensor, part of the space environment monitor system aboard GOES, consists of 
two ion chamber detectors, which provide whole-sun X-ray fluxes in the 0.05-0.3 and 0.1-0.8 
nm wavelength bands. Solar soft X-ray flares are classified according to their peak burst 
intensity measured in the 0.1-0.8 nm wavelength band by GOES. The letters (A, B, C, M, 
X) denote the order of magnitude of the peak flux on a logarithmic scale, and the number 
following the letter gives the multiplicative factor, i.e.. An = n x 10~^, Bn = n x 10~^, 
Cn = n X 10~^, Mn = n x 10"^ and Xra = n x 10~^ W/m^. In general, n is given as a 
float number with one decimal (prior to 1980, n is listed as an integer). No background 
subtraction is applied to the data. Now the data is widely available from the NOAA Space 



Environment Center site (http://www.ngdc.noaa.gov/stp/SOLAR/ftpsolarflares.html). 



In the meantime, a wealth of data has been accumulated. It makes worthwhile re- 
investigating the tempora l and s patial features of soft X-ray (SXR) flares on an extensive 



statistical basis. iLi et al.l (119981 ) studied the distribution of the X-ray flares (M>1) from 
1987 to 1992 with respect to helio longitude. They have shown that the flares were not 
uniformly distributed in longitude. The temporal analy sis of X-flare statis t ics concerns 



basically the waiting-time distribution (see, for example, iBoffetta et al.l Il999l : iMoon et al. 
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200ll : iLepreti et al. boOll : IWheatland booi IVeronig et al.l |2002| and so on). In the present 



analysis we make use of SXR flares observed by GOES during 1976-2006. Our consideration 
will be devoted only to the energy statistics of soft X-ray solar flares in time. 



4. Predictive Tool Description 

Our analysis is based on the propertie s of fractional autoregressive integrated moving 



average (FARIMA) processes (iBeran Ill994l ). They are widely used in modeling of various 
complex physical systems. The FARlMA{p,d,q) process is defined as the solution of the 
equation ^{B)A'^X{n) = 6(i?)e„, n E Z, where B is the shift operator BX{n) = X{n — 1) 
and A is the difference operator, i. e. AX{n) = X{n) — X{n — 1). Here $ and G are 
the polynomials of degree p and q respectively, d takes fractional values, either positive or 
negative, and "innovations" ej are independent and identically distributed (i.i.d.) random 
variables. The polynomials $ and correspond to autoregressive (AR) and moving average 
(MA) parts, respectively. The linear representation of FARIMA processes takes the form 

oo 

^in) = ^Cn-jej , (1) 



for details see (iBeran 111994 ). The innovations may be either Gaussian, non-Gaussian with 
finite variance or they may have infinite variance. For infinite variance innovations e, one 
may consider, for example, symmetric and skewed stable distributions, as well as Pareto 
distributions. Both are characterized by the parameter a and their tails P{e > x) satisfy 

P{e > x) = 1 - F{x) r-^ x'"', as a; oo, (2) 

where F{x) denotes the corresponding distribution function and ~ denotes that the ratio of 
the left-hand side to the right-hand one tends to 1, as x — oo. It should be noted that the 
Levy-stable distributions have < a < 2 whereas for the Pareto distribution the parameter 
a is greater than zero. The resulting process X{n) will be long-range dependent and Levy- 
stable if the innovations are Levy-stable, and asymptotically will be in the domain of attrac 



tion of a Levy-stable distribution if the innovations are Pareto (see lSamorodnitsky fc Taqqu 



19941 ). Moreover, such FARIMA processes are asymptotically self-similar with d — 1/a. 



The Levy-stable distribution, named after the French mathematician Paul Levy who 
investigated the behavior of sums of independent random variables, is most conveniently 
described by its characteristic function (f){6) - the inverse Fourier transform of the probability 
density function. The most popular form of the characteristic function of a Levy-stable 
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random variable is given by the expression 



log 



—a 



-i/?sign(^)tan(7ra/2)} + 
a\9\{l + 2ip sign(^) log \9\/7i} + i^O, 



a ^ 1, 
1, 



a 



(3) 



wher e < a < 2, — 1 < /3 < 1, a > and yU G R are parameters of this distribu- 



tion (ISamorodnitsky fc Taqqu Ill994l ). The Pareto distribution, introduced by the Italian 
economist Vilfredo Pareto, is a power law probability density tha t we represent in the form 
of f{x) = aX"{X + a;)~"~^, where A and a are positive constants (IBurnecki et al.ll2005l ). 



The power-law behavior of the tails implies that the variance is infinite if a < 2. The 
tail index (exponent) a controls the rate of decay of the tail of the distribution function F. 
Modeling with FARIMA time series with infinite variance allows to take into account heavy 
tails. Through a suitable choice of coefficients Cn-j one can also add long-term memory ef- 
fects. The FARIMA processes is an useful family of models because it offers a lot of flexibility 
in modeling long-range and short-range dependence by choosing the memory parameter d 
and appropriate autoregressive and moving average coefficients in expression ([1]). 

The problem of estimating the exponent in heavy-tailed data has a long history in 
statistics because of its p ractical imp ortance. The presence of heavy tails in data was 
firstly noted in the work of IZipf 1 ( 19321 ) in his study of word frequencies in languages. Next, 
Mandelbrot I(jl960l ) noted their presence in financial data. Since the early 1970s the heavy- 
t ailed behavior has been noted in rnany o ther scientific fields (see, for example, reviews of 
Adler et al.lll998l : iPark fc Willinger 1120001 ). However, the availability of huge amount of var- 
ious data poses a set of new challenges for the problem of estimating the tail index. The 
point is that the data can be contaminated by extraneous oscillations, different noises with 
f inite variance and so o r i. This makes the analysis of heavy-tailed data more complicated 
( jjanicki fc Weron I1994I : iLynch et al.l 120051) . The t ime series of soft X-ray solar emission 
relates to such problematic data (iBaiesi et al. Therefore, for reliability we will esti- 

mate the tail index by different statistical tests. One of them is based on the a symptotic 
max self- similarity properties of heavy-tailed maxima (jStoev fc Michailidis In this 

test the maximum values of data are calculated over blocks of size m, scaled at rate of m^/". 
By examining a sequence of growing block sizes m = 2^ , 1 < j < loggA^, j G N, and sub- 
sequently estimating the mean of logarithms of block-maxima one obtains an estimation of 
the tail index a. Another estima tor, that we use, under the assumption of the Levy stable 
law applies the iMcCuUoch I(ll986l ) quantile fit. 
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5. Cycling of Self-similarity 

Using the max self-similarity estimator and considering our data in year intervals, we 
have analyzed how the solar cycling influences on the tail index. Figure [1] shows a clear 
correlation between solar activity and the index. When the solar activity is around maxima, 
the tail index is larger than one, whereas in minima it tends to fall down less t han one. The 



index value in the period of high solar activity almost coincides with the result of lWeron et al. 



( I2OO5I ). Although their data insert only X-ray solar bursts of C and M type, the value a was 
estimated to be 1.2674. The present analysis extends the tail index analysis on some cycles 
and speaks surely that the index tendency observed earlier is kept at least during the recent 
three solar cycles. 

The McCuUoch's testing dMcCulloch Il986h gives similar results. This test shows that 



the tail parameter a depends on the solar activity value during the three solar cycles. Of 
course, the estimator has a particular character. It is convenient for the analysis of random 
variables with a stable distribution. Nevertheless, the distribution also has heavy tails. Our 
flrst aim is to flnd such distribution that will be best fltted to the experimental series. 



6. Solar flares 

Now we restrict our attention to such time intervals, in which the solar activity is 
strong (in particular, 1978-1981, 1988-1992 and 1999-2003). We use X-ray flare data from 
GOES satellite, that contain information about time of appearance and energy of solar flares 
(from http://www.ngdc.noaa.gov/stp/SOLAR/ftpsolar-flares.html). The captured energy 
was transmitted by X-rays emitted during blasts on a solar surface from 2000 January 1 to 
2002 December 31. We aggregated the energy values on a daily basis. The time series is 
presented in Fig. [2J 

The flrst estimation procedure of the self-similarity parameter H (i. e. the Hurst ex- 
ponent)is the so-called flnite impulse response transformation (FIRT). The FIRT estimator 
involve an array of coefficients. The array is made out of flnite impulse response coeflicients. 
The estimator Hpi^x is obtained by per forming a log-linear regression on the coefficients and 
measuring the slope (jStoev et al. It is important to note that the estimator Hf 



I FIRT 



is unbiased for all a falling in the range (0, 2). 



An alternative method of testing sc aling and correla tion properties of a time series is 



the variance of residuals method (VR) (jPeng et al.lll994j ). First, the series is divided into 



blocks of size m. Then, within each block, the partial sums of the series are calculated. A 
least-squares line is fitted to the partial sums within each block, and the sample variance of 
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the residuals is computed. The variance of residuals is proportional to m^^ . This variance 
of residuals is computed for each block, and the median (or average) is computed over the 
blocks. A log-log plot versus m should follow a straight line with a slope of 2H. 

The R/S method is one of t he oldest and better known methods for estimation of 
the Hurst parameter (IHurst Il95lh . Hurst found that drought in the Nile Valley is not a 
random phenomenon, but rather that the region is inclined to become progressively more 
arid after a succession of long droughts. It i s widely known that for Ga ussian (i.e. finite 
variance) time series the method returns H (IMandelbrot fc Wallis Il969h . The popularity 
of the method has been also a source of misunderstandings and errors. This is due to the 
fact that, in general, for power-law distributed time series (i.e. with infinite variance) the 
method yields d + 1/2. In particular for Pareto or Levy stable distributions the output is 
H — 1/ a + 1/2 (iTaqqu &: Teverovsky Il998h . The method is based on R/S (rescaled adjusted 
range) statistics. The series is divided into blocks. Then, within each block, the statistics 
is calculated. Finally, arithmetic means of the values of the statistics over the blocks are 
calculated and a least squares line is fitted to the mean for different lengths of the blocks. 
The slope should be equal to d+ 1/2. 

For the finite variance cases, the interpretation of the FIRT and VR estimators is very 
similar to the Hurst exponent: if only short-range correlations (or no correlations at all) 
exist in the studied series, then Hfirt = Hvr = 1/2; if there is a correlation then Hfifit = 
HvR 7^ 1/2. Moreover, if the estimator Hfirt = Hvr is greater than 1/2, the time series is 
persistent and if Hfirt = Hyu < 1/2, then the time series is not persistent. 

Note that both estimators give an information on memory and not on distribution of 
the process increments. 

The analysis of the data shows that the tails of the underlying distribution conform to 
the power law. Hence, we model the data by a FARIMA process with Pareto innovations. 
As the power-law dist ributions belong to the domain of attraction of stable law (see, eg. 
Janicki &: Weron Il994h the resulting distribution of the FARIMA process should be close 
to the stable one. We applie d the McCuUoch quantile fit to obtain the parameters of the 
distribution (IMcCulloch Il986l ) . The value of a was estimated to be 1.213, see Figured One 
may check that the estimated value of a of simulated FARIMA times series with Pareto 
innovations is usually underestimated, see e. g. Figure [31 Therefore, we assume that the 
innovations in our model follow the Pareto law with a = 1.25. 



According to I Weron et al.l (120051 ). in order to recover both the self-similarity exponent H 
and the memory parameter d (hence, the distribution parameter a) we can use the following 
BMW^ computer test. The surrogate data are obtained here by a random shuffling of the 
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original data positions. 



• If the process is FARIMA with Gaussian noise, then the values of the estimator should 
change to 1/2 for the surrogate data independently on the initial values. 

• If the process is FARIMA with a-stable or Pareto noise for a < 2, then the values of 
the estimator should change to 1/a for the surrogate data independently on the initial 
values. 



The above formalism can be easily applied to determine basic features of an empirical 
data series. Now we employ this to study an empirical time series recorded from the system 
describing the energy of solar flares (Fig. [2]). The obtained values of the parameters are listed 
in Table [TJ Therefore, from the results for the surrogate data, the corresponding estimates 
for the parameter 1/a are: l/apiRT = Hfirt = 0.8452 and l/avR = Hvr = 0.7722. 
We observe that the estimators are close to the one assumed in our model: 1/aMC = 0.8. 
Moreover, we choose ci = 0.19 as the highest admissible value of d, which is close to the one 
obtained via the R/S method for the original data, see Tabled! 

One may notice that the estimators of H obtained via FIRT and VR methods, see 
Table [1], are great er than theoret i cally a dmissible in the FARIMA model, i. e. they exceed 



one. As stated in iBurnecki et al.l (l2008l ). this can be justified by performing simulations of 
the FARIMA processes and estimating the parameter H on the simulated time series via 
different methods. It appears that a reasonable percentage of the values of H is higher than 
the data estimates, so one can not reject the hypothesis that the underlying model is the 
FARIMA(0,(i,0) process. Nevertheless, we now look for an enhanced model which would de- 
scribe better the behavior of different estimators obtained for the original and shuffled data, 
and would improve the fit in terms of the prediction error. Thus, we propose a slight gen- 
eralization of the model incorporating the short- dependence component, namely FARIMA 
{2,d,0) model. We estimated the AR(2) coefficients: ai (linear term) and 02 (quadratic 
term) via the mean-square error (MSE) minimalization scheme taking into account three 
statistics: FIRT, VR and RS, see Figure HI The estimated values are: ai = 0.02 and 
02 = 0.03. The FARIM A processes were generated according to the algorithm presented by 



Stoev fc Taqqu \t004 \ 



We calculate Hfirt, HyR and djis estimators for the simulated FARIMA {2,d,0) pro- 
cesses (top panel in Figs. O and [H]) and the corresponding shuffied data (bottom panel in 
Figs. Oandin]). We generate 1000 trajectories of size 2^° which is close to the length of the 
original solar data, i. e. 1089 and present the results in form of the so-called box plots. The 
box plot has lines at the lower quartile, median, and upper quartile values. The whiskers 
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are lines extending from each end of the box to show the extent of the rest of the data. 
Outhers are data with values beyond the ends of the whiskers. If there is no data outside 
the whisker, a dot is placed at the bottom whisker. The results correspond to the analysis 
of the original solar data included in Table [TJ Thus, our FARIMA simulations reconstruct 
well the structure of the original solar flares data. 

We conclude that the proper model could be based on the FARIMA(2,(i,0) process with 
Pareto innovations with the parameters Oi = 0.02, 02 = 0.03, d = 0.19, and a = 1.25 



which has the long-range dependence property since d > 1 — 2/ a (IBurnecki et al.ll2008l ). The 
Pareto distribution is also convenient for the present analysis because it gives a description 
of positive random variables whereas the Levy-stable one for 1 < a < 2 is related to both 
positive and negative random variables, but a time series of x-ray flare energy is quite positive 
by definition. 

We also calculated the 1-day-ahead prediction for the FARIMA f2,(j,0) time ser ies (for 
the prediction discussion in the infinite variance FARIMA case see iKokoszka Il995[ l. The 
results are depicted in Figure [71 



7. Concluding Remarks and Discussion 

In this paper we demonstrate how self-similar models driven by Levy stable noise 
can be useful for modeling X-ray solar data. To be more precise we have suggested the 
FARIMA(2,(i,0) model with a-stable noise for predicting solar flare appearance in the pe- 
riod of a strong solar activity. 

The procedure is illustrated in Section [6] for the captured energy transmitted by X- 
rays emitted during blasts on a solar surface from 2000 January 1 to 2002 December 31. 
Comparing the values of the different estimators for the original data series and for the 
surrogate data we estimated the components of the self-similarity index corresponding to the 
memory of the time series (d) and to the tail properties of the time series values distribution 
(a), see Tabled! Thus, this allows in principle to build a proper physical model for analyzing 
the solar activity. 

The analysis of soft X-ray emission observations shows that this series is enough com- 
plicated in nature. We have seen a strong dependence of statistics on solar cycling in the 
data. However, if one takes year intervals, then this cycle influence becomes not so strong. 
This allows one to reconstruct a statistical model for predicting the soft X-ray solar activity 
on the nearest solar cycle. The soft X-ray solar time series contains both long-range de- 
pendence and heavy-tailed effects. The first creates a random number of strong bursts on 
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a background, and the second forms their persistence between each other. The most conve- 
nient model for their joint description is a FARIMA time series. In view of solar cycling the 
model permits one to predict a time series of soft X-ray solar flares, when the solar activity 
will be again near its maximum in 2010-2014 years. While wc do not claim that this model 
provides the only possible explanation, it does provide a rigorous statistical picture of the 
expected observations of X -ray solar flares. 

A.A.S. is grateful to the Institute of Physics and the Hugo Steinhaus Center for Stochas- 
tic Methods for pleasant hospitality during his visit in Wroclaw University of Technology. 
The GOES X-ray hght curve was made available courtesy of the NOAA Space Environment 
Center, Boulder, CO. 
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Fig. 1. — The evolution of the tail index a during the last solar cycles 1974-2006 (top 
picture) obtained via the max self-similarity method, the bottom picture shows Wolf numbers 
(characterizing solar activity) in this period. 
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Fig. 2. — Energy-time series of solar flares from 2000 January 01 to 2002 December 31. 
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Fig. 3. — Values of the calculated a estimators for the simulated FARIMA time series. Solid 
line represents the value of the estimator for the analyzed data. 
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Fig. 4. — Mean squared error on the basis of the calculated H and d estimators for the 
simulated FARIMA time series with respects to the estimators obtained for the analyzed 
data for different ai and 02- 
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Fig. 5. — Values of the FIRT and VR estimators for the simulated time series (top panel) 
and the surrogate data (bottom panel) for the generated FARIMA processes. Solid line 
represents the value of the estimator for the analyzed data. 
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Fig. 6. — Values of the RS estimator Hjis for the simulated time series (top panel) and the 
surrogate data (bottom panel) for the generated FARIMA processes. Solid line represents 
the value of the estimator for the analyzed data. 
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Fig. 7. — Solar flare data and prediction in the FARIMA(2,d,0) model. 
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Table 1: Values of the FIRT, VR and RS estimators for the original time series and the 
shuffled (surrogate) solar flare data. 



Data set Hfirt Hvr 


dRS 


Original time series 


Solar flares 1.1424 1.0665 


0.2408 


Surrogate data 


Solar flares 0.8452 0.7722 


0.0507 



